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Abstract 

We calculate the ratio of the viscosity to the entropy density for both Bose and Fermi 
gases in the unitary limit using a new approach to the quantum statistical mechanics 
of gases based on the S-matrix. In the unitary limit the scattering length diverges and 
the S-matrix equals —1. For the fermion case we obtain 77/5 > 4.7 times the proposed 
lower bound of h/AirkB which came from the AdS/CFT correspondence for gauge theories, 
consistent with the most recent experiments. For the bosonic case we present evidence that 
the gas undergoes a phase transition to a strongly interacting Bose-Einstein condensate, 
and is a more perfect fluid, with rj/s > 1.3 times the bound. 



I. INTRODUCTION 



One measure of the perfection of a fluid is its viscosity. In recent years new 
insights on such properties came from string theory, more precisely the AdS/CFT 
correspondence [l, 2|. This correspondence relates a conformally invariant (scale in- 
variant) strongly coupled gauge theory in + 1 space-time dimensions to a gravita- 
tional dual in one higher spatial dimension. By studying black hole solutions in the 
higher dimensional theory, one can study finite temperature and density properties 
of the lower dimensional quantum field theory. In its original version, the conformal 
quantum field theory is a certain = 4 supersymmetric gauge theory, which is con- 
formally invariant for all couplings since the beta function vanishes. Although such 
a supersymmetric gauge theory does not describe nature as we currently understand 
it, the AdS/CFT correspondence is nevertheless very useful for thinking about these 
difficult problems in new ways and inspiring the study of specific properties that were 
hardly considered before. A prominent example is the ratio of the shear viscosity t] 
to the entropy density s. In natural units with h = = 1, it is a. dimensionless 
quantity. Using the AdS / CFT correspondence it was found that the supersymmetric 
gauge theory had rj/s = 1/Att. It was conjectured that this value represents a lower 
bound ^, i.e. 

n > LA. (1) 

s An Kb 

The observed rj/s for the quark-gluon plasma, studied via heavy ion collisions, is 
about 5 times this bound jsj. Theoretical estimates for the quark-gluon plasma tend 
to be lower, but consistent with the bound. 

Note that the bound ([1]) does not depend on the speed of light. Thus, the natural 
question arises: "Do non-relativistic condensed matter systems respect the conjec- 
tured bound, and which class of fiuids has the smallest 77/5?" For the relativistic 
fiuids studied wth the AdS/CFT correspondence, the scale invariance plays an im- 



portant role. 



n 



i has thus been proposed that the unitary Fermi gas may be the most 



perfect fluid j#t7|, in part because it is scale invariant. In the so-called unitary limit 
of a quantum Bose or Fermi gas, the scattering length a diverges. This occurs at 
a fixed point of the renormalization group, thus these systems provide interesting 
examples of interacting, scale-invariant theories with dynamical exponent 2; = 2, i.e. 
non-relativistic. They can be realized experimentally bytuning the scattering length 
to ±00 using a Feshbach resonance. (See for instance jsl, l9| and references therein.) 
They are also thought to occur at the surface of neutron stars. These systems have 
attracted much theoretical interest 1CH24|. 

Because of the scale-invariance, the only length scales in the problem are based on 



the density n^^^, and the thermal wavelength = ^J2'K /mT. Equivalently, the only 
energy scales are the chemical potential /i and the temperature T. The problem is 
challenging since there is no small parameter to expand in such as no?. Any possible 
critical point must occur at a specific value of x = fi/T. This can be translated into 
universal values for ncAji, or for fermions universal values for T^/Tp where ep = ksTp 
is the Fermi energy. For instance the critical point of an ideal Bose gas is the simplest 
example, where ricX^ = C(3/2) = 2.61. 

The models considered are the simplest models of non-relativistic bosons or 
fermions with quartic interactions. The bosonic model is defined by the action for a 
complex scalar field 0. 



s 



d xdt 



2m 



(2) 



For fermions, due to the fermionic statistics, one needs at least a 2-component field 



a=t,i 



(3) 



In both cases, positive g corresponds to repulsive interactions. The bosonic theory 



only has a f/(l) symmetry. The fermionic theory on the other hand has the much 
larger SO (5) symmetry. This is evident from the workj24| which considered an N- 
component version with Sp(2N) symmetry, and noting that Sp(4) = S0(5). 

The original AdS/CFT conjecture was for a specific, relativistic, and supersym- 



metric gauge theory. There have been some proposal 
spondence to learn about non-relativistic systems |25l- 



s to use the AdS/CFT corre- 
. One difficulty is that the 



conformal symmetry of relativistic systems is larger than the Schrodinger symme- 
try of non-relativistic systems, so the black hole solutions on the gravity side are 
less obvious. Also, given a black hole geometry, it remains unclear which condensed 
matter system it is dual to. Thus far, the AdS/CFT approaches for non-relativistic 
systems lead to rj/s = 1/An. Recent experimental work on the unitary Fermi gas 



reports values oi rj/s about 4 — 5 times the bound(32j, thus it seems unlikely that a 
gravity dual exists that corresponds exactly to the unitary Fermi gas. Nevertheless, 
it is hoped that one can still discover some general, model-independent properties, in 
the same way that AdS/CFT for supersymmetric Yang-Mills provided insights into 
QCD. 

For the remainder of this paper we will describe a novel, but more conventional 



approach to the problem [291]. The main approximation made is that we only consider 
2-body interactions, and consistently resum their contributions to the free energy 
via an integral equation. The 2-body interactions are expressed in terms of the zero- 
temperature S-matrix, which can be calculated exactly. The viscosity is calculated 
in a kinetic, semi-classical approach. The calculation predicts a minimum rj/s oi 
4.7 times the conjectured lower bound. This is consistent with the most recent 



experiments j32j, which suggests our approximation is a good one. 

Theoretical studies have mainly focussed on the fermionic case, and for the most 
part at zero temperature, which is appropriate for a large Fermi energy. The bosonic 
case has been less studied, since a homogeneous bosonic gas with attractive interac- 



tions is thought to be unstable against mechanical collapse, and the collapse occurs 
before any kind of BEC. The situation is actually different for harmonically trapped 



gases, where BEC can occur SSj. However studies of the homogeneous bosonic case 



were based on a small, negative scattering length 34N37l|. and it is not clear that the 



conclusions reached there can be extrapolated to the unitary limit. Since the den- 



sity of collapse is proportional to l/a[35|, extrapolation to infinite scattering length 
suggests that the gas collapses at zero density, which seems unphysical, since the gas 
could in principle be stabilized at finite temperature by thermal pressure. One can 
also point out that in the van der Waals gas, the collapse is stabilized by a finite size 
of the atoms, which renders the compressibility finite. In the unitary limit, there is 
nothing to play such a role. In the sequel we will present evidence that the unitary 
Bose gas undergoes BEC when riA^ ~ 1.3. This lower value compared to the free 
case is consistent with the attractive interactions. For this bosonic case, the mini- 



mum rj/s predicted in 29| was only 1.3 times the bound, and thus may be a better 



candidate than unitary fermions for the most perfect strongly interacting fluid. 

II. S-MATRIX, RENORMALIZATION GROUP AND SCATTERING 
LENGTH 

In this section we describe the renormalization group fixed point, the bound state, 
and the scattering length. The interplay between all of these properties is most clearly 
seen from the S-matrix. These results are well-known and can be found in the book 



30| for instance. We include them here for completeness and in order to consistently 
describe all of our conventions. 

The free versions of the models described in the introduction have a scale invari- 
ance with dynamical exponent z = 2, i.e. are invariant under 



As we now explain, the models possess a renormalization group fixed point, i.e. 
quantum critical point, where they have the same scale invariance. The renormal- 
ization group behavior can be inferred from the coupling constant dependence of the 
S-matrix. Consider first the single boson; the differences for two-component fermions 
will be described at the end of this section. The S-matrix can be calculated exactly 
by summing multi-loop ladder diagrams. By the Galilean invariane, the two-body S- 
matrix depends only on the difference of the incoming momentum of the two particles 
k,k': 

^' 167i/mgji + i\k-k'\ ^ ' 

Unitarity of the S-matrix amounts to S*S = 1. 

The momentum space integrals for the higher loop corrections are divergent and an 
upper cut-off A must be introduced. In the above expression, gji is the renormalized 
coupling: 

1 1 mA 

— = - + (5) 
9r 9 47r^ 

Defining g = ^/A, where ^ is dimensionless, and requiring gji to be independent of 
A gives the beta-function: 

where i = —log A is the logarithm of a length scale. The above beta function is 
exact since it was calculated from the exact S-matrix. One thus sees that the theory 
possesses a fixed point at the negative coupling g^, = —An'^/mK. 

We turn now to the scattering length a. From the above expression for the S- 
matrix one can infer the scattering amplitude and compute the total cross-section a. 
Equating a = no? gives 

a{k) = ^ , (7) 
27r ^1 + {mgRk/^Ti f 
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where here k is the momentum of one of the particles in the center of mass frame. 
If a{k) is measured at very small momentum transfer |k — k'| ^0, this leads to the 
definition of the scattering length 

mgu, mg 

a = = -, ; — r 8 

2n 271(1 ~g/g,) ^' 

One sees that scattering length diverges at precisely the fixed point g = g^. Note 

the S-matrix becomes S* = — 1. The scattering length a — ±oo, depending on from 

which side g* is approached. When g = g~ , i.e. just less than g^, then a — )■ oo, 

whereas when g = g^ , a ^ —oo. 

Finally, we turn to the bound state. The S-matrix (j4]) has a pole at k = IQ-Ki/mgR. 

Since physical bound states correspond to poles at Im(A;) > 0, the bound state exists 

only for g below The energy of this bound state is 

1287r2 

-Abound— state 5 o~ y^j 

Note that at the fixed point where g^ diverges, the energy of the bound state goes 
to zero as it should, since it disappears beyond this point. 

Consider now two-component fermions with the action in the Introduction. The 
relative normalizations of the coupling g were chosen such that the beta function 
is the same for both the boson verses fermion cases. The S-matrix eq. @ here 
represents the scattering of two fermions of opposite spin. Thus, the fermion case 
also has a diverging scattering length at the fixed point. In the fermionic context, the 
coupling g^ is the boundary of the so-called BEC/BCS crossover. On the BCS side 
just above g^, the scattering length is negative, which implies effectively attractive 
interactions. Here it is believed that at low enough temperatures there is a phase 
transition to a strongly interacting version of superconductivity. For g < g^,, the 
scattering length is positive, signifying repulsive interactions, and a bound state 
exists. This bosonic bound state can undergo Bose-Einstein condensation, hence 



this region is referred to as the BEC side. The physics is expected to be smooth 
as one crosses (7*. In the treatment of the thermodynamics below, we will work 
on the BCS side since here there is no need to incorporate a bound state into the 
thermodynamics. 



III. THERMODYNAMICS AT THE QUANTUM CRITICAL POINT 

At the quantum critical point, the only energy scales in the problem are the 
chemical potential fi and the temperature T = 1//3. This implies some universal 



131]. The free energy density 



scaling forms for the various thermodynamic functions 
has the form 

jr = -C(5/2)T A^-^ c(^/T) (10) 



where ( is Riemann's zeta function and = \/2'k jmT is the thermal wavelength. 
The scaling function c is only a function of x = /i/T. With the above normalization, 
a single free boson has c = 1 in the limit of x — > 0. It is also convenient to define the 
scaling function g, which is a measure of the quantum degeneracy, in terms of the 
density n as follows: 

nA^ = q (11) 

The two scaling functions c and q are of course related since n = —dT/dfi, which 
leads to g = ({5/2)c', where c' is the derivative of c with respect to x. Henceforth b' 
will always denote the derivative of a function b{x) with respect to x. 

The approach to the statistical mechanics of particles developed in comple te g en 
is based on the S-matrix. On starts from a formula derived in 



erality in 



38 



for 



the partition function: 

Z = Zo + — [ dEe-^^ Tr Im^e log S(E) (12) 
Airi J 
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where Zq is the partition function for the free theory and S{E) is the off-shell S- 
matrix operator in the usual scattering theory. Though the above expression is 
simple enough, a considerable amount of additional work is required to turn it into 
something useful. The trace is over the multi-particle Fock space, thus the above 
expression contains contributions from N-body processes for all N. For integrable 
theories in one spatial dimension, the N-body S-matrix factorizes into 2-body S- 
matrices, and though it has never been proven, the above expression should recover 
the thermodynamic Bethe ansatz 40 1 . 
As described in detail in 138 



one can develop a diagrammatic description of the 
contributions to eq. f[T^ . where vertices with 2N legs correspond to the logarithm of 
the S-matrix for N particle scattering and lines connection the vertices are occupation 
numbers. The free energy can be obtained from a variational principle based on 
a Legendre transformation between the chemical potential and the density. The 
variational principle leads to an integral equation satisfied by the occupation numbers 
with kernels involving the logarithm of the S-matrix. In the present context, we 
restrict ourselves to 2-body processes only, which, as explained in the last section, 
can be calculated exactly. This should be a good approximation if the gas is not too 
dense. We also remark that this method is quite different than the method used for 
example by Nozieres-Schmitt-Rink although both involve the S, or t-matrix. The 
latter formalism does not lead to the self-consistent integral equations we describe 
below, but rather a BCS-like equation. We point out that this latter method was 
extended to the unitary gas in \3l\ . 



Let us now describe the main result derived in [38| . Consider again for simplicity 
a single component bosonic or fermionic gas. The filling fractions, or occupation 
numbers, are parameterized in terms of a pseudo-energy e(k): 



9 



which determine the density: 



/ 1 .... 

^ / (27r)3 e/^^'M -s ^ ' 



where s = 1, — 1 corresponds to bosons, fermions respectively. The consistent sum- 
mation of 2-body scattering leads to an integral equation for the pseudo-energy e(k), 
analogous to the Yang- Yang integral equation. It is convenient to define the quantity: 

where = k^/2m. Then y satisfies the integral equation 

^(k) = l + /3y (^^(k-kOjy^^ (16) 



The free energy density is then 



J" 



(27r)3 



-s logfl — se 



i(i-r 



17) 



2 e^' 

The kernel G is related to the logarithm of the 2-body S-matrix of the last section, 
and depends on the coupling g. In the unitary limit g ^ g* the kernel simplifies 
greatly since S = —1: 

where the — sign corresponds to g being just below the fixed point g^,, where the 
scattering length a — )■ +oo on the BEC side, whereas the -|- sign corresponds to 
a — 7- — oo on the BCS side. As explained in the last section, we will work on the BCS 
side. 

Finally comparing with the definitions above for the scaling functions c, q one 
finds: 

q{x) = ^ r dnV^^^y^^^ (19) 
Jo e« - sy{K)z 
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and 

^ ^ [ ""^ - - ^^^^) 

where z = e^l^ is the fugacity and the dimensionless integration variable is k = 
k^/2mT. The ideal, free gas limit corresponds to y = 1 where q = sLi^/2{sz) and 
c = sLi^/2{sz)/({5/2), where Li is the polylogarithm. The BEC critical point of the 
ideal gas occurs at /i = 0, i.e. q = C(3/2). 

Consider now two-component fermions with the action ([3]). There are two pseudo- 
energies e^^i satisfying two coupled integral equations. Due to the SU(2) symmetry, 
for equal chemical potentials = e^. However the available phase space for 2-particle 
scattering is doubled and the kernels have an extra 1/2: 

G^fcrmi — '^Gyjose (21) 

Thus the two-component fermion reduces to two identical copies of the above 1- 
component expressions, with the modification ( I2T]) . 



IV. PHASE TRANSITIONS IN TEMPERATURE 

In the present context of quantum statistical mechanics, the quantum critical 
point at g = does not signify a phase transition. However in the unitary limit at 
fixed g = g^f phase transitions may occur in regions of temperature versus density. 
From the scaling form of the free energy, such a phase transition must occur at a fixed, 
specific value of x = Xc = yu/T. We will continue to refer to such phase transitions as 
critical points. The simplest case is BEC in the free boson theory, where the phase 
transition occurs at Xc = 0, which can be translated into the well-known relation 
between temperature and density at the critical point rZcAj, = C(3/2) = 2.61. 



11 



A. Fermionic case 



The integral equation for eq. ( IT6|) . can be solved numerically by iteration. 
One first substitutes = 1 on the right hand side and this gives the approximation 
Hi for y. One then substitutes yi on the right hand side to generate y2, etc. For 
regions of z where there are no critical points, this procedure converges rapidly, and 
as little as 5 iterations are needed. For fermions, as one approaches zero temperature, 
i.e. X large and positive, more iterations are needed for convergence. The following 
results are based on 50 iterations. 

When z <^ 1, y ^ 1, and the properties of the free ideal gas are recovered, since the 
gas is very dilute. There are solutions to eq. flTB]) for all — oo < x < oo. {x = ^/T). 
The scaling function c at zero chemical potential is c(0) = 1.76 compared to the free 
fermion value Cfree = 2 — 2~^/^ = 1.29, thus the interactions have a significant effect. 
(These are twice the 1-component values.) 

Whereas c and q are nearly featureless, other quantities seem to indicate a phase 
transition at large density. For instance, the entropy per particle decreases with de- 
creasing temperature up to x < ~ 11.2. Beyond this point the entropy per particle 
has the unphysical behavior of increasing with temperature. A further indication that 
the region x > Xc is unphysical is that the specific heat per particle becomes nega- 
tive, as shown in Figure [H When x ^ 0, CyjN approaches the classical value 3/2. 
This leads us to suggest a phase transition, a.t x = Xc = 11.2. This value of x can be 
expressed as a critical temperature Tc in units of the Fermi energy ep = ksTp- From 
the definition ep = (3-7r^n/-\/2)^^'^/'n2, one has T/Tf = (4/3v/7rg)2/3. From the value 
q{xc), one finds the critical temperature T^/Tp ^ 0.1. As we will show, our analysis 
of the viscosity to entropy-density ratio suggests a higher Tc/Tp- There have been 
numerous estimates of Tc/Tp based on various approximation schemes, mainly using 



Monte Carlo methods on the lattice 



17 



22| . quoting results for T^/Tp between 0.05 



12 



and 0.23. The work [18| puts an upper bound T^Tp < 0.14, and the most recent 
results of Burovski et. al. quote T^/Tp = 0.152(7). Our result is thus consistent 
with previous work. 



Cv/N 



Specific heat per particle for fermions 




X = fi/T 



FIG. 1: Specific heat per particle as a function of x for fermions. 



B. Bosonic case 

The possibility of a phase transition in the unitary Bose gas is more subtle, since 
a bosonic gas with attractive interactions is susceptible to mechanical collapse to a 
denser state. This issue has been studied for small negative scattering length in a 



number of works 



34 



-37l|. and the consensus is that mechanical collapse occurs before 
BEG. However it is not clear whether this conclusion can be extrapolated to unitary 
limit where the scattering length is infinite. In the approach described in the last 
section, we found strong evidence for a phase transition to a strongly interacting 
version of BEG, as we now explain. 
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We again solved the integral equation (fT6!) by iteration, starting from y = 1. 
Since the occupation numbers decay quickly as a function of n, we introduced a cut- 
off k < 10. For X less than approximately —2, the gas behaves nearly classically. The 
main feature of the solution to the integral equation is that for x > Xc = —1.2741, 
there is no solution that is smoothly connected to the classical limit x — t- — oo. 
Numerically, when there is no solution the iterative procedure fails to converge. In 
Figure [2], we plot £:(k = 0) as a function of x, and one sees that it goes to zero at 
Xc- This implies the occupation number / diverges at k = at this critical point. 
One clearly sees this behavior in Figure [3l We also found that the compressibility 
diverges at Xc, again consistent with BEG. 



e(k = 0)/T 




^ X = jj/T 



FIG. 2: The pseudo-energy e at k = as a function of x = fi/T. 

This strongly suggests that there is a critical point at Xc which is a strongly 
interacting, scale invariant version of the ideal BEG. In terms of the density, the 
critical point is: 

ncA^ = 1.325, (/i/T = Xc = -1.2741) (22) 
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FIG. 3: The occupation number /(k) for x = —1.275 and Xc = —1.2741. 

The negative value of the chemical potential is consistent with the effectively attrac- 
tive interactions. The above should be compared with the ideal BEC of the free 
theory, where Xc = and rZcA^ = C(3/2) = 2.61, which is higher by a factor of 2. 
A critical exponent v characterizing the diverging compressibility can be defined 

as 

n^{T-T,Y^ (23) 

A log-log plot of the compressibility verses T — Tc shows an approximately straight 
line, and we obtain u ^ 0.69. This should be compared with BEC in an ideal gas, 
where ~ 1.0. Clearly the unitary gas version of BEC is in a different universality 
class. 
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V. VISCOSITY TO ENTROPY DENSITY RATIO 



Finally we turn to the intended focus of this article, the ratio of the viscosity to 
entropy density. Consider first a single component gas. The simplest expressions for 
the shear viscosity are based on kinetic theory, where it is related to the momentum 
transfer through an imaginary 2-dimensional plane cutting through the 3-dimensional 
bulk in the presence of a velocity gradient in the fluid flow. It can be expressed as 

7] = ^nvmifT-ee (24) 



where v is the average speed and ^frcc is the mean free pathj41|. The mean free path 
is £free = 1/ (V^na) where a is the total cross-section. (The \/2 comes from the ratio 
of the mean speed to the mean relative speed.) 

Consider flrst the bosonic model. (There were two errors by factors of 2 in the 



original discussion presented in 29| which are here corrected. However they turn out 



to cancel and the main result below eq. (1251) is the same as before. ) In the unitary 



limit where the S-matrix S = —1, the amplitude is 



29|: 



The cross section was calculated by the standard manner in which it is related to an 
integration over flnal state momenta and the square of the amplitude Ai 

Stt |k-k'|2 ^ ' 

Finally, 

(|k-kf) = (k^ + k'2 - 2k ■ k') =2mV (27) 

V ' ' / avc V / avc ^ ' 

since the cosine of the angle between k and k' averages to zero. This gives 

m V 
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(30) 



For a unitary gas, the relation between the energy, volume and pressure is the 
same as for a free theory: E/V = 3p/2^^. Since the pressure is due to the kinetic 
energy, this implies 

lmv' = E/N = ^^T (29) 

Thus 7] oc {cT/c')^^'^. In the high temperature limit, the function y — )■ 1 and c oc 

Li5/2{e^) oc X and c' oc Li3/2(e^) oc x as x — )■ 0. Thus, as expected r] oc T'^/^. Since 

the entropy density s = —dJ^/dT, one finally has 

r] _ VStt /c\3/2 1 
s ~ 8C(5/2) \d) 5c/2 - xd 

For two-component fermions, the available phase space, referred to as X in 29 1 
is doubled. Also, spin up particles only scatter with spin down. This implies the 
cross section is 8 times smaller than the expression (I26p . Since the entropy density is 
doubled due to the two components, this implies that 77/5 is 4 times the expression 
eq. (1301) . where c is the 1-component value appropriate to fermions. 

The ratio rj/s for fermions as a function of T/Tp is shown in Figure HJ and is 
in good agreement both quantitatively and qualitatively with the experimental data 
summarized in 61]. The lowest value occurs at x = 2.33, which corresponds to 
T/Tp = 0.28, and 

^ > 4.72-4- (31) 

s AirkB 

This is consistent with the most recent experimental data which shows a minimum 



that is about 4 — 5 times the bound 



32|. Other, considerably more complicated 



theoretical approaches, e.g. using a Kubo formula for the viscosity, gives results 
around 6 — 7 times the bound 42 1. 



For bosons, the ration 77/5 is plotted in Figure [5] as a function of T/T^. One sees 
that it has a minimum at the critical point, where 

a > 1.26-^ (32) 

s AixkB 
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Viscosity/entropy for attractive fermions 
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FIG. 4: The viscosity to entropy-density ratio as a function of T/Tp for fermions. The 
horizontal line is I/Att. 

Thus the bosonic gas at the unitary critical point is a more perfect fluid than that 
of fermions. 

VI. CONCLUSIONS 



In summary, we have calculated the ratio of viscosity to entropy density for both 
Bose and Fermi gases in the unitary limit using the S-matrix based approach in 29 1. 
For the fermionic gas, we found rj/s > is greater than 4.7 times the conjectured 
lower bound of l/47r, which is lower than the result of other theoretical approaches. 



and more consistent with the most recent experimental results[32j. We provided 
evidence that the unitary Bose gas is stable and has a strongly interacting BEC 
phase transition. The same calculation indicates that the Bose case is a more perfect 
fluid, with rj/s > 1.7 times the lower bound. 
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Viscosity to entropy ratio for bosons 




FIG. 5: The viscosity to entropy-density ratio as a function of T/Tc for bosons. Tlie 
iiorizontal line is l/47r. 
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